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FOREWORD 


Methodology  and  Characterization  studies  during  fiscal  years 
1977  and  1978  (FY  77  and  78)  included  gravity  surveys  in  ten 
valleys  in  Arizona  (five),  Nevada  (two),  New  Mexico  (two),  and 
California  (one).  The  gravity  data  were  obtained  for  the  pur¬ 
pose  of  estimating  the  gross  structure  and  shape  of  the  basins 
and  the  thickness  of  the  valley  fill.  There  was  also  the  possi¬ 
bility  of  detecting  shallow  rock  in  areas  between  boring  loca¬ 
tions.  Generalized  interpretations  from  these  surveys  were  in¬ 
cluded  in  Fugro  National's  Characterization  Reports  (FN-TR-26a 
through  e) . 

During  the  FY  77  surveys,  measurements  were  made  to  form  an 
approximate  1-mile  grid  over  the  study  areas  and  contour  maps 
showing  interpreted  depth  to  bedrock  were  made.  In  FY  79,  the 
decision  was  made  to  concentrate  on  verifying  and  refining 
suitable  area  boundaries.  This  decision  resulted  in  a  reduction 
in  the  gravity  program.  Instead  of  obtaining  gravity  data  on  a 
grid,  the  reduced  program  consisted  of  obtaining  gravity 
measurements  along  profiles  across  the  valleys  where  Verifica¬ 
tion  studies  were  also  performed. 

The  Defense  Mapping  Agency  (DMA),  St.  Louis,  was  requested 
to  provide  gravity  data  from  their  library  to  supplement  the 
gravity  profiles.  For  Big  Smoky,  Hot  Creek,  and  Big  Sand 
Springs  valleys,  a  sufficient  density  of  library  data  are  avail¬ 
able  to  permit  construction  of  interpreted  contour  maps  instead 
of  just  two-dimensional  cross  sections. 

In  late  summer  of  FY  79,  supplementary  funds  became  available  to 
begin  data  reduction.  At  that  time,  inner  zone  terrain  correc¬ 
tions  were  begun  on  the  library  data  and  the  profiles  from  Big 
Smoky  Valley,  Nevada,  and  Butler  and  La  Posa  valleys,  Arizona. 
The  profile  data  from  Whirlwind,  Hamlin,  Snake  East,  White 
River,  Garden,  and  Coal  valleys,  Nevada,  became  available  from 
the  field  in  early  October  1979. 

A  continuation  of  gravity  interpretations  has  been  incorporated 
into  the  FY  80  program,  and  the  results  are  being  summarized  in 
a  series  of  valley  reports.  Reports  covering  Nevada-Utah 
gravity  studies  will  be  numbered  "FN-TR-33-"  followed  by  the 
abbreviation  for  the  subject  valley.  In  addition,  more  detailed 
reports  of  the  results  of  FY  77  surveys  in  Dry  Lake  and  Ralston 
valleys,  Nevada,  are  being  prepared.  Verification  studies  are 
continuing  in  FY  80,  and  gravity  studies  are  included  in  the 
program.  DMA  will  continue  to  obtain  the  field  measurements, 
and  it  is  planned  to  return  to  the  grid  pattern.  The  interpre¬ 
tation  of  the  grid  data  will  allow  the  production  of  contour 
maps  which  will  be  valuable  in  the  deep  basin  structural  analy¬ 
sis  needed  for  computer  modeling  in  the  water  resources  program. 

i 


FN-TR-33-BS 


The  gravity  interpretations  will  also  be  useful  in  Nuclear  Hard¬ 
ness  and  Survivability  (NH&S)  evaluations. 

The  basic  decisions  governing  the  gravity  program  are  made  by 
BMO  following  consultation  with  TRW,  Inc.,  Fugro  National,  and 
the  DMA.  Conduct  of  the  gravity  studies  is  a  joint  effort 
between  DMA  and  Fugro  National.  The  field  work,  including 
planning,  logistics,  surveying,  and  meter  operation  is  done  by 
the  Defense  Mapping  Agency  Hydrographic/Topographic  Center 
(DMAHTC) ,  headquartered  in  Cheyenne,  Wyoming.  DMAHTC  reduces 
the  data  to  Simple  Bouguer  Anomaly  (see  Section  A1 . 4 ,  Appendix 
A1.0).  The  Defense  Mapping  Agency  Aerospace  Center  (DMAAC) ,  St. 
Louis,  Missouri,  calculates  outer  zone  terrain  corrections. 

Fugro  National  provides  DMA  with  schedules  showing  the  valleys 
with  the  highest  priorities.  Fugro  National  also  recommended 
locations  for  the  profiles  in  the  FY  79  studies  with  the  pro¬ 
vision  that  they  should  follow  existing  roads  or  trails.  Any 
required  inner  zone  terrain  corrections  are  calculated  by  Fugro 
National  prior  to  making  geologic  interpretations. 
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1.0  INTRODUCTION 


1 . 1  OBJECTIVE 

Gravity  measurements  were  made  in  Big  Smoky  Valley  for  the 
purpose  of  estimating  the  overall  shape  of  the  structural  basin, 
the  thickness  of  alluvial  fill,  and  the  location  of  concealed 
faults.  The  estimates  will  be  useful  in  modeling  the  dynamic 
response  of  ground  motion  in  the  basin  and  in  evaluating  ground- 
water  resources. 

1.2  LOCATION 

Big  Smoky  Valley  is  in  northeastern  Esmeralda  and  northwestern 
Nye  counties,  Nevada.  The  town  of  Tonopah,  Nevada,  (Figure  1) 
is  near  the  southeast  corner  of  the  valley  on  highways  6  and  95. 
Access  throughout  the  valley  is  good  due  to  an  extensive  network 
of  well  maintained,  unpaved  roads.  The  valley  area  is  princi¬ 
pally  used  for  rangeland,  but  effects  of  extensive  mining  opera¬ 
tions  are  apparent  around  the  valley  edge. 

Big  Smoky  Valley  (Figure  2)  is  bounded  on  the  east  by  t he  San 
Antonio  mountains  and  the  Toquima  Range,  on  the  west  by  Rjyston 
Hills  and  Monte  Cristo  Range,  on  the  south  by  Lone  Mountain  and 
Tonopah,  Nevada,  and  on  the  north  by  the  Toiyabe  Range.  The 
area  covered  by  this  report  lies  betweer  latit  .  ies  J8'!uo'  and 
38045’  and  longitudes  117o00l  and  117o40l.  It  is  approximately 
40  miles  by  51  miles  (64  km  by  81  km). 
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1.3  SCOPE  OF  WORK 

A  total  of  1386  gravity  stations  were  used  in  this  report.  The 
Defense  Mapping  Agency  Aerospace  Center  ( DMAAC )  supplied  1157 
gravity  stations  from  their  library  and  229  new  gravity  sta¬ 
tion  measurements  were  made  by  the  Defense  Mapping  Agency 
Hydrographic- Topographic  Center/Geodetic  Survey  Squadron 
( DMAHTC/GSS ) . 
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2.0  GRAVITY  Dr»TA  REDUCTION 


DMAHTC/GSS  obtained  the  basic  observations  for  the  new  stations 
and  reduced  them  to  Simple  Bouguer  Anomalies  (SBA)  as  described 
in  Appendix  A1.0.  Up  to  three  levels  of  terrain  corrections 
were  applied  to  the  new  stations  to  convert  the  SBA  to  the 
Complete  Bouguer  Anomaly  (CBA) .  Only  the  first  two  levels  of 
terrain  corrections  described  below  were  applied  to  the  library 
stations . 


First,  the  DMAAC ,  St. Louis,  Missouri,  used  its  library  of 
digitized  terrain  data  and  a  computer  program  to  calculate 
corrections  out  to  104  miles  (167  km)  from  each  station.  When 
the  program  could  not  calculate  the  terrain  effects  near  a 
station,  a  ring  template  was  used  to  estimate  the  effect  of 
terrain  within  approximately  3000  feet  (914  m)  of  the  station. 
The  third  level  of  terrain  corrections  was  applied  to  those 
stations  where  10  feet  (3  m)  or  more  of  relief  was  observed 
within  130  feet  (40  m)  .  In  these  cases,  the  elevation  differ¬ 
ences  were  measured  in  the  field  at  a  distance  of  130  feet 
along  six  directions  from  the  stations.  These  data  were 
used  to  calculate  the  effect  of  the  very  near  relief. 
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3.0  GEOLOGIC  SUMMARY 


Big  Smoky  Valley  and  the  general  vicinity  lies  in  an  area  con¬ 
sidered  to  be  within  the  upper  plate  of  the  Roberts  Mountain 
thrust  sheet  along  the  western  limit  of  the  miogeosyncl inal 
belt.  The  combination  of  this  older  al lochthonous  sequence  with 
Tertiary  volcanic  rocks  and  tectonic  events  result  in  complex 
geologic  structure  and  extremely  diverse  rock  types  in  the 
mountains  surrounding  Big  Smoky  Valley.  To  the  east,  the  west¬ 
ern  flank  of  the  San  Antonio  Mountains  is  composed  of  Tertiary 
volcanic  flows  (rhyolites,  andesites,  and  basalts)  and  tuffs 
with  some  sediments,  predominantly  of  the  Esmeralda  Formation 
(sandstones,  siltstones,  and  mudstones)  (Kleinhampl  and  Ziony, 
1967).  The  southern  Toquima  Range,  at  the  northeastern  boundary 
of  Big  Smoky  Valley,  consists  of  Cretaceous  (granites)  and 
Tertiary  rhyolite,  basalt,  and  tuffs  with  some  metamorphic  and 
sedimentary  rock  (Ferguson  and  Cathcart,  1954;  Kleinhampl  and 
Ziony ,  1967) . 

The  Toiyabe  Range  at  the  northern  end  of  the  site  is  chiefly 
Tertiary  air-fall  and  ash-flow  tuff  with  minor  outcrops  of 
Paleozoic  shale  and  limestone  along  the  eastern  flank.  Rocks 
along  the  western  side  of  Big  Smoky  Valley  (Royston  Hills  and 
Monte  Cristo  Range)  are  chiefly  Tertiary  tuffs,  rhyolites, 
andesites,  and  basalts. 

Basin-fill  deposits  within  the  valley  reach  combined  thicknesses 
of  about  3000  feet  (915  m)  (Rush  and  Schroer,  1970).  Younger 
alluvial  fans,  commonly  in  association  with  fluvial,  eolian,  and 
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playa  deposits,  occur  throughout  the  central  valley,  locally 
extending  to  the  mountain  fronts.  Deposits  consist  predominant¬ 
ly  of  sand  with  increasing  gravel  content  near  steep  mountain 
fronts,  I n t e rmed i a t e- 1  eve  1  alluvial  fans  are  restricted  to 
mountain  front  areas  and  are  rarely  exposed  along  the  valley 
axis.  Stream  channel  and  playa  deposits  occur  in  northern  and 
central  Big  Smoky  Valley  along  major  south  flowing  drainages. 

Prominant  north  to  northeast  trending  range-bounding  faults 
occur  on  the  steep  eastern  side  of  the  Toiyabe  Range  and  on  the 
western  side  of  the  San  Antonio  Mountains.  Evidence  of  late 
Quaternary  faulting  is  seen  in  the  form  of  fault  scarps  in  allu¬ 
vium  along  the  southwestern  edge  of  Lone  Mountain,  along  the 
entire  eastern  side  of  the  valley,  and  at  scattered  locali¬ 
ties  throughout  the  valley.  The  faults  exhibiting  greatest 
apparent  displacement  occur  along  the  west  side  of  the  San 
Antonio  Mountains  where  scarps  are  75  feet  (23  m)  or  more  in 
height.  Many  faults  previously  mapped  along  the  eastern  side 
(Rush  and  Schroer,  1979)  do  not  show  displacement  of  the  allu¬ 
vium  but  may  be  represented  as  tonal  and  vegetation  lineaments. 
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4.0  INTERPRETATION 

The  basis  of  interpretation  is  the  Complete  Bouguer  Anomaly 
(CBA)  shown  in  Drawing  1.  The  CBA  is  defined  in  Appendix 
A1.4. 

The  CBA  data  are  reduced  to  a  set  of  values  at  the  points  of  a 
uni formly-spaced  geographic  array,  or  grid.  This  facilitates 
the  mathematical  treatment  of  irregularly  spaced  data.  Gravity 
stations  are  numerous  along  roads  and  areas  where  access  is 
easy.  The  grid-process  is  done  using  an  algorithm  which  com¬ 
putes  a  value  at  each  grid  point  from  the  gravity  station  data 
within  a  circular  area  around  the  grid  point.  A  bell-shaped 
weighting  function  assigns  greater  weight  to  the  nearer  data 
points.  The  grid-point  spacing  is  chosen  to  match  the  average 
data  spacing.  A  1.2  mile  (2  km)  grid  spacing  was  used  for  this 
analysis . 

4.1  REGIONAL-RESIDUAL  SEPARATION 

A  fundamental  step  in  gravity  interpretation  is  isolation  of 
the  part  of  the  CBA  which  represents  the  geologic  feature  of 
interest,  in  this  case  the  valley  fill.  The  valley  fill  has  a 
lower  density  than  the  bedrock  and,  therefore,  creates  a  nega¬ 
tive  gravity  anomaly.  The  portion  of  the  CBA  which  corresponds 
to  this  alluvial  material  is  called  the  "residual  anomaly.” 

The  CBA  contains  long-wave-length  components  from  deep  and  broad 
geologic  structures  extending  far  beyond  the  valley.  These 
long-wave-length  components,  called  the  regional  gravity,  have 
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been  approximated  by  a  second-degree  trend  surface.  The  trend 
surface  is  a  least-squares  fit  to  the  CBA  values  at  the  bedrock 
stations  of  the  region.  This  regional  trend  is  subtracted  from 
the  CBA  to  obtain  the  residual  anomaly.  The  residual  anomaly  is 
used  to  calculate  a  simple  geologic  model  which  fits  the  gravity 
data  and  is  consistent  with  geologic  knowledge  from  other 
sources . 

4.2  DENSITY  SELECTION 

The  construction  of  a  geologic  model  from  the  residual  anomaly 
requires  selection  of  density  values  representative  of  the  allu¬ 
vial  fill  and  of  the  underlying  rock.  Since  only  very  general¬ 
ized  density  information  is  available,  the  geologic  interpreta¬ 
tion  of  the  gravity  data  can  only  be  a  coarse  approximation. 
Average  in  situ  density  of  the  fill  material  was  measured  in  six 
shallow  borings  (Table  1)  at  a  depth  of  161  feet  (49  m)  .  The 
observed  density  range  for  the  soil  was  1.7  to  2.3  g/cm^.  The 
largest  measured  density  value  was  used  in  the  modeling  process, 
instead  of  the  average,  because  the  overall  alluvium  density  is 
expected  to  increase  due  to  compaction  with  depth  (compaction 
with  depth  and  age  is  discussed  by  Woollard,  1962  and  Grant  and 
West,  1965) . 

The  basement  material  underlying  the  Big  Smoky  basin  is  thought 
to  be  the  Paleozoic  carbonate  rocks  which  are  found  sparsely  in 
the  surrounding  mountain  ranges.  Published  values  for  carbonate 
rocks  typically  range  between  2.6  and  2.8  g/cm^.  The  Paleozoic 
carbonate  rocks  in  Nevada  are  generally  reported  to  be  relative¬ 
ly  high  in  density,  on  the  order  of  2.8  g/cm3#  This  value  was 
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selected  to  represent  the  density  of  the  basement  rock*  The 
density  contrast  used  for  modeling  was  -0.50  g/cm^. 

4.3  MODELING 

Modeling  was  done  with  the  aid  of  a  computer  program  which 
iteratively  calculates  a  three-dimensional  solution  of  gravity 
anomaly  data  (Cordell,  1970).  The  gravity  anomaly  is  represent¬ 
ed  by  discrete  values  on  a  two-dimensional  grid.  The  source  of 
the  anomaly  (the  volume  of  low-density  valley  fill)  is  repre¬ 
sented  by  a  set  of  vertical  prism  elements.  The  tops  of  the 
prisms  lie  in  a  common  horizontal  plane.  The  bottoms  of  the 
prisms  collectively  represent  the  bottom  of  the  valley  fill. 
Each  prism  has  a  cross-sectional  area  equal  to  one  grid  square 
and  a  uniform  density.  A  grid  square  of  1.2  miles  by  1.2  miles 
(2  km  by  2  km)  was  selected  as  representative  of  the  gravity 
station  distribution.  Computations  were  made  for  five  itera¬ 
tions  of  mutually  interactive  prism  adjustment.  The  root-mean- 
square  error  for  the  entire  grid  was  less  than  0.5  milligal. 

The  calculated  thickness  of  the  valley  fill  depends  upon  the 
accuracy  of  the  -0.50  g/cm^  density  contrast  (i.e.,  fill  den¬ 
sity  minus  rock  density)  used.  Since  neither  density  is  per¬ 
fectly  known,  nor  even  uniform,  the  calculated  thickness  should 
be  expected  to  contain  a  corresponding  degree  of  uncertainty. 
Table  1  lists  three  seismic  refraction  lines  and  two  deeper 
borings  which  were  used  as  constraints  in  the  modeling  process. 
The  seismic  refraction  lines  located  near  the  mountain  flanks 
recorded  high  velocities  which  may  represent  the  basement 
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material  . 

The  alluvial 

fill 

material 

in  the  center 

o  f 

the 

valley  is 

approx imately 

1400 

to  1500 

feet  (427  to 

457 

m) 

thick  according  to  the  two  literature  borings  (Table  1).  The 
calculated  thickness  of  fill  or  interpreted  depth  to  rock  is 
contoured  in  Drawing  2. 

4.4  DISCUSSION  OF  RESULTS 

The  depth-to-rock  contours  (Drawing  2)  depict  a  long  north-south 
graben  with  multiple  normal  faults.  The  interpretation  was 
aided  by  regional  geologic  information  from  published  reports, 
aerial  photographs,  and  field  reconnaissance  which  located  the 
surface  expression  of  faults.  The  second  vertical  derivative  of 
the  CB A  field  was  calculated  to  guide  the  placement  of  the 
faults  shown.  Because  the  zero  contour  of  the  second  vertical 
derivative  marks  the  steepest  part  of  the  input  CBA  field,  the 
faults  were  placed  along  the  zero  contour. 

The  fault  system  on  the  west  side  of  the  San  Antonio  Mountains 
and  Toquima  Range  is  generally  close  to  the  mountain  front.  The 
fault  system  is  also  close  to  the  mountain  front  on  the  east 
side  of  the  Toiyabe  Range.  However,  there  appears  to  be  a  pedi¬ 
ment  several  miles  wide  east  of  Royston  Hills  and  Monte  Cristo 
Range.  At  the  southern  end  where  the  valley  bifurcates,  faults 
are  interpreted  to  be  near  both  flanks  of  Lone  Mountain. 

Several  northwest-southwest  trending  fault  zones  are  transverse 
to  the  major  faults.  There  is  no  indication  in  regional  or  lo¬ 
cal  geologic  data  or  in  this  gravity  data  that  these  features 
are  related  to  through-going  tectonic  systems.  They  appear  to 


CRO 


IO! 


IC 


FN-TR-33-BS 


13 


represent  zones  of  differential  movement  within  the  central 
down-dropped  block. 

The  northwest  trending  grabens  between  the  Toiyabe  Range  and  the 
Royston  Hills  and  between  the  Toquima  Range  and  San  Antonio 
Mountains,  however,  may  represent  fundamental  crustal  fracture 
systems  but  their  role  in  the  regional  tectonic  framework  is  not 
understood  at  this  time  and  is  beyond  the  scope  of  this  report. 

In  summary,  central  Big  Smoky  Valley  appears  to  comprise  a  junc¬ 
tion  of  several,  variously  oriented  horsts  and  grabens.  The 
central  area  is  a  small  graben  which  is  distinct  from  grabens  in 
northeastern  Big  Smoky  Valley  and  those  to  the  south  on  both 
sides  of  Lone  Mountain.  This  northeast-trending  major  graben 
system  is  disrupted  by  northwest-trending  minor  grabens. 
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5.0  CONCLUSIONS 


The  interpretation  of  Big  Smoky  Valley  gravity  data  indicates 
that  there  are  major  range-bounding  normal  fault  systems  extend¬ 
ing  the  full  length  of  the  valley  on  both  sides.  The  maximum 
depth  calculated  in  the  valley  is  4000  feet  (1219  m)  . 

The  calculated  bedrock  depths  are  only  approximations  because 
little  is  known  about  the  actual  density  distribution  in  and 
around  the  valley  and  the  residual  gravity  anomaly  is  necessar¬ 
ily  based  on  an  interpreted  regional  field.  An  average  density 

/  • 

contrast  of  0.50  g/cm^  between  the  alluvium  and  bedrock  was 
used  to  calculate  the  thickness  of  the  valley-fill  material. 
Future  studies  that  acquire  better  density  data  or  measure 
actual  depths  to  bedrock  in  deep  parts  of  the  valley  can  be  used 
to  refine  the  gravity  interprepation . 


-fuGRO  inc 


FN-TR-33-BS 


15 


BIBLIOGRAPHY 


Cordell,  Lindreth,  1970,  Iterative  solution  of  three-dimensional 
gravity  anomaly  data,  Geol.  Survey  Computer  Contr.  No.  10: 
U.S.  Geol.  Survey,  Washington,  D.C. 

Ferguson,  H.  G.  and  Cathcart,  S.  H.  (1954)  Geology  of  the  Round 
Mountain  quadrangle,  Nevada:  U.  S.  Geol.  Survey,  Geol. 

Quad.  Map  GQ-40 ,  scale  1:125,000. 

Fugro  National,  Inc.  1979,  MX  siting  Investigation  geotechnical 
evaluation,  Verification  studies,  v.  IB,  FN-TR-27. 

,  1979,  Verification  study  data  of  Big  Smoky  CDP, 

Nevada,  v.  VIII,  FN-TR-27. 

Kleinhampl,  F.  J.,  and  Ziony,  J.  I.  1967,  Preliminary  geologic 
map  of  northern  Nye  County,  Nevada:  U.  S.  Geol.  Survey, 

Open  file  map,  scale  1:200,000. 

Rush,  F.  E.,  and  Schroer,  C.  V.,  1979,  Water  resources  of  Big 

Smoky  Valley,  Landu,  Nye  and  Esmeralda  counties,  Nevada: 
Nev.  Dept,  of  Conserv.  and  Nat.  Resources  Bull.  No.  41, 
84  p. 

Stewart,  John  H.  1980,  Geology  of  Nevada:  Nevada  Bureau  of 
Mines  and  Geology,  Special  Pub.  4. 

Stewart,  John  H.  and  Carlson,  John  E.  1978,  Geology  of  Nevada: 
Nevada  Bureau  of  Mines  and  Geology. 


IGRO  N 


1C 


FN-TR-3  3-BS 


APPENDIX  A 1 . 0 

GENERAL  PRINCIPLES  OF  THE 
GRAVITY  EXPLORATION  METHOD 


IGRO  NATIONAL  INC 


FN-TR-33-BS 


A1.0  GENERAL  PRINCIPLES  OF  THE  GRAVITY 
EXPLORATION  METHOD 

Al.J  GENERAL 

A  gravity  survey  involves  rnea s u r emen t  of  differences  in  the 
g r a v  i  t a t i o na 1  field  between  various  points  on  the  earth's 
surface.  The  gravitational  field  values  being  measured  are  the 
same  as  those  influencing  all  objects  on  the  surface  of  the 
earth.  They  are  generally  associated  with  the  force  which 
causes  a  1  gm  mass  to  be  accelerated  at  980  cm/sec2.  This 
force  is  normally  referred  to  as  a  1-g  force. 

Even  though  in  many  applications  the  gravitational  field  at  the 
earth's  surface  is  assumed  to  be  constant,  small  but  distin¬ 
guishable  differences  in  gravity  occur  from  point  to  point. 
In  a  gravity  survey,  the  variations  are  measured  in  terms 
of  milligals.  A  milligal  is  equal  to  0.001  cm/second2  or 
0.00000102  g.  The  differences  in  gravity  are  caused  by  geo¬ 
metrical  effects,  such  as  differences  in  elevation  and  latitude, 
and  by  lateral  variations  in  density  within  the  earth.  The 
lateral  density  variations  are  a  result  of  changes  in  geologic 
conditions.  For  measurements  at  the  surface  of  the  earth,  the 
largest  factor  influencing  the  pull  of  gravity  is  the  density  of 
all  materials  between  the  center  of  the  earth  and  the  point  of 
measurement . 

To  detect  changes  produced  by  differing  geological  conditions, 
it  is  necessary  to  detect  differences  in  the  gravi tat ional  field 
as  small  as  a  few  milligals.  To  recognize  changes  due  to 
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geological  conditions,  the  measurements  are  "corrected"  to  ac¬ 
count  for  changes  due  to  differences  in  elevation  and  latitude. 

Given  this  background,  the  basic  concept  of  the  gravitational 
exploration  method,  the  anomaly,  can  be  introduced.  If,  instead 
of  being  an  oblate  spheroid  character i zed  by  complex  density 
variations,  the  earth  were  made  up  of  concentric,  homogeneous 
shells,  the  gravitational  field  would  be  the  same  at  all  points 
on  the  surface  of  the  earth.  The  complexities  in  the  earth's 
shape  and  material  distribution  are  the  reason  that  the  pull  of 
gravity  is  not  the  same  from  place  to  place.  A  difference  in 
gravity  between  two  points  which  is  not  caused  by  the  effects  of 
known  geometrical  differences,  such  as  in  elevation,  latitude, 
and  surrounding  terrain,  is  referred  to  as  an  "anomaly." 

An  anomaly  reflects  lateral  differences  in  material  densities. 
The  gravi ta tionai  attraction  is  smaller  at  a  place  underlain  by 
relatively  low  density  material  than  it  is  at  a  place  underlain 
by  a  relatively  high  density  material.  The  term  "negative 
gravity  anomaly"  describes  a  situation  in  which  the  puli  of 
gravity  within  a  prescribed  area  is  small  compared  to  the  area 
surrounding  it.  Low-density  alluvial  deposits  in  basins  such  as 
those  in  the  Nevada-Utah  region  produce  negative  gravity  anoma¬ 
lies  in  relation  to  the  gravity  values  in  the  surrounding 
mountains  which  are  formed  by  more  dense  rocks. 

The  objective  of  gravity  exploration  is  to  deduce  the  variations 
in  geologic  conditions  that  produce  the  gravity  anomalies 
identified  during  a  gravity  survey. 
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A  1.2  INSTRUMENTS 

The  sensing  element  of  a  LaCoste  and  Romberg  qravimeter  is  a 
mass  suspended  by  a  zero-length  spring.  Deflections  of  the 
mass  from  a  null  position  are  proportional  to  changes  in  gravi¬ 
tational  attraction.  These  instruments  are  sealed  and  compen¬ 
sated  for  atmospheric  pressure  changes.  They  are  maintained  at 
a  constant  temperature  by  an  internal  heater  element  and  thermo¬ 
stat.  The  absolute  value  of  gravity  is  not  measured  directly  by 
a  gravimeter.  It  measures  relative  values  of  gravity  between 
one  point  and  the  next.  Gravitational  differences  as  small  as 
0.01  milligal  can  be  measured. 

A 1 . 3  FIELD  PROCEDURES 

The  gravimeter  readings  were  calibrated  in  terms  of  absolute 
gravity  by  taking  readings  twice  daily  at  nearby  USGS  ora^itv 
base  stations.  Gravimeter  readings  fluctuate  becau^  of  small 
time-related  deviations  due  to  the  effect  of  earth  tides  and 
instrument  drift.  Field  readings  were  corrected  to  account  for 
these  deviations.  The  magnitude  of  the  tidal  correction  was 
calculated  using  an  equation  suggested  by  Goguel  (1954): 

C  =  P  +  Ncos  &  (cos  0  +  sin  0  )  +  Scos  ^  (cos  -  sin  0) 
where  C  is  the  tidal  correction  factor,  P,  N,  and  S  are  time- 
related  variables,  and  0  is  the  latitude  of  the  observation 
point.  Tables  giving  the  values  of  P,  N,  and  S  are  published 
annually  by  the  European  Association  of  Exploration  Geophysi¬ 
cists  . 
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The  meter  drift  correction  was  based  on  readings  taken  at  a 
designated  base  station  at  the  start  and  end  of  each  day.  Any 
difference  between  these  two  readings  after  they  were  corrected 
for  tidal  effects  was  considered  to  have  been  the  result  of 
instrumental  drift.  It  was  assumed  that  this  drift  occurred  at 
a  uniform  rate  between  the  two  readings.  Corrections  for  drift 
were  typically  only  a  few  hundredths  of  a  milligal.  Readings 
corrected  for  tidal  effects  and  instrumental  drift  represented 
the  observed  gravity  at  each  station.  The  observed  gravity 
values  represent  the  total  gravitational  pull  of  the  entire 
earth  at  the  measurement  stations. 

A1.4  DATA  REDUCTION 

Several  corrections  or  reductions  are  made  to  the  observed 
gravity  to  isolate  the  portion  of  the  gravi  tat ional  pull  which 
is  due  to  the  crustal  and  near-surface  materials.  The  gravity 
remaining  after  these  reductions  is  called  the  "Bouguer  Anoma¬ 
ly."  Bouguer  Anomaly  values  are  the  basis  for  geologic  inter¬ 
pretation.  To  obtain  the  Bouguer  Anomaly,  the  observed  gravity 
is  adjusted  to  the  value  it  would  have  had  if  it  had  been 
measured  at  the  geoid,  a  theoretically  defined  surface  which 
approximates  the  surface  of  mean  sea  level.  The  difference 
between  the  "adjusted”  observed  gravity  and  the  gravity  at  the 
geoid  calculated  for  a  theoretically  homogeneous  earth  is 
the  Bouguer  Anomaly. 
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Four  separate  reductions,  to  account  for  four  geometrical 
effects,  are  made  to  the  observed  gravity  at  each  station  to 
arrive  at  its  Bouguer  Anomaly  value. 

a.  Free-Air  Effect;  Gravitational  attraction  varies  inversely 
as  the  square  of  the  distance  from  the  center  of  the  earth. 
Thus,  corrections  must  be  applied  for  elevation.  Observed  grav¬ 
ity  levels  are  corrected  for  elevation  using  the  normal  vertical 
gradient  of: 

FA  =  -0.09406  mg/ft  (-0.3086  mi 1 1 igals/meter ) 
where  FA  is  the  free-air  effect  (the  rate  of  change  of  gravity 
with  distance  from  the  center  of  the  earth).  The  free-air 
correction  is  positive  in  sign  since  the  correction  is  opposite 
the  effect. 

b .  Bouguer  Effect :  Like  the  free-air  effect,  the  Bouguer 

effect  is  a  function  of  the  elevation  of  the  station,  but  it 
considers  the  influence  of  a  slab  of  earth  materials  between 
the  observation  point  on  the  surface  of  the  earth  and  the 
corresponding  point  on  the  geoid  (sea  level).  Normal  practice, 
which  is  to  assume  that  the  density  of  the  slab  is  2.67  grams 
per  cubic  centimeter  was  followed  in  these  studies.  The  Bouguer 
correction  (Bc),  which  is  opposite  in  sign  to  the  free-air 
correction,  was  defined  according  to  the  following  formula. 

Bc  =  0.01276  (2.67)  hf  (milligals  per  foot) 

Bc  =  0.04185  (2.67)  hm  (milligals  per  meter) 
where  hf  is  the  height  above  sea  level  in  feet  and  hm  is  the 
height  in  meters. 
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c.  Latitude  Effect:  Points  at  different  latitudes  will  have 
different  "gravities”  for  two  reasons*  The  earth  (and  the 
geoid)  is  spheroidal,  or  flattened  at  the  poles.  Since  points 
at  higher  latitudes  are  closer  to  the  center  of  the  earth  than 
points  near  the  equator,  the  gravity  at  the  higher  latitudes  is 
larger.  As  the  earth  spins,  the  centrifugal  acceleration 
causes  a  slight  decrease  in  gravity.  At  the  higher  latitudes 
where  the  earth's  radii  are  smaller,  the  centrifugal  accel¬ 
eration  diminishes.  The  gravity  formula  for  the  Geodetic 
Reference  System,  1967,  gives  the  theoretical  value  of  gravity 
at  the  geoid  as  a  function  of  latitude.  It  is: 

g  =  978.0381  (1  +  0.00532 04  sin2  ff  -  0.0000058  sin2?jz*)  gals 
where  g  is  the  theoretical  acceleration  of  gravity  and  0  is 
the  latitude  in  degrees.  The  positive  term  accounts  for  the 
spheroidal  shape  of  the  earth.  The  negative  term  adjusts  for 
the  centrifugal  acceleration. 

The  previous  two  corrections  (free  air  and  Bouguer)  have  ad¬ 
justed  the  observed  gravity  to  the  value  it  would  have  had  at 
the  geoid  (sea  level).  The  theoretical  value  at  the  geoid  for 
the  latitude  of  the  station  is  then  subtracted  from  the  adjusted 
observed  gravity.  The  remainder  is  called  the  Simple  Bouguer 
Anomaly  (SBA) .  Most  of  this  gravity  represents  the  effect  of 
material  beneath  the  station,  but  part  of  it  may  be  due  to 
irregularities  in  terrain  (upper  part  of  the  Bouguer  slab)  away 
from  the  station. 
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d.  Terrain  Effect:  Topographic  relief  around  the  station  has 
a  negative  effect  on  the  gravitational  force  at  the  station.  A 
nearby  hill  has  upward  gravitational  pull  and  a  nearby  valley 
contributes  less  downward  attraction  than  a  nearby  material 
would  have.  Therefore,  the  corrections  are  always  positive. 
Corrections  are  made  to  the  SBA  when  the  terrain  effects  were 
0.1  milligal  or  larger.  Terrain  corrected  Bouguer  values  are 
called  the  Complete  Bouguer  Anomaly  ( CB A ) .  When  the  CB A  is 
obtained,  the  reduction  of  gravity  at  individual  measurement 
points  (stations)  is  complete. 

A 1 . 5  INTERPRETATION 

The  first  step  in  interpretation  is  to  separate  the  portion  of 
the  CBA  that  might  be  caused  by  the  light-weight,  basin-fill 
material  overlying  the  heavier  bedrock  material  which  forms  the 
surrounding  mountains  and  presumably  the  basin  floor.  Since  the 
valley-fill  sediments  are  absent  at  the  stations  read  in  the 
mountains,  the  CBA  values  at  these  bedrock  stations  are  used  as 
the  basis  for  constructing  a  regional  field  over  the  valley.  A 
regional  field  is  an  estimation  of  the  values  the  CBA  would  have 
had  if  the  light-weight  sediments  (the  anomaly)  had  not  been 
there  . 

The  difference  between  the  CBA  and  the  regional  field  is 
called  the  "residual”  field  or  residual  anomaly.  The  residual 
field  is  the  interpreter's  estimation  of  the  gravitational 
effect  of  the  geologic  anomaly.  The  zero  value  of  the  residual 
anomaly  is  not  exactly  at  the  rock  outcrop  line  but  at  some 
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distance  on  the  "rock"  side  of  the  contact.  The  reason  for  this 
is  found  in  the  explanation  of  the  terrain  effect.  There  is  a 
component  of  g ravi tat ional  attraction  from  material  which  is 
not  directly  beneath  a  point. 

If  the  "regional”  is  well  chosen,  the  magnitude  of  the  residual 
anomaly  is  a  function  of  the  thickness  of  the  anomalous  (fill) 
material  and  the  density  contrast.  The  density  contrast  is  the 
difference  in  density  between  the  alluvial  and  bedrock  material. 
If  this  contrast  were  known,  an  accurate  calculation  of  the 
thickness  could  be  made.  In  most  cases,  the  densities  are  not 
well  known  and  they  also  vary  within  the  study  area.  In  these 
cases,  it  is  necessary  to  use  typical  densities  for  materials 
similar  to  those  in  the  study  area. 

If  the  selected  average  density  contrast  is  smaller  than  the 
actual  density  contrast,  the  computed  depth  to  bedrock  wi  1 1  be 
greater  than  the  actual  depth  and  vice-versa.  The  computed 
depth  is  inversely  proportional  to  the  density  contrast.  A  ten 
percent  error  in  density  contrast  produces  a  ten  percent  error 
in  computed  depth.  An  iterative  computer  program  is  used  to 
calculate  a  subsurface  model  which  will  yield  a  gravitational 
field  to  match  (approximately)  the  residual  gravity  anomaly. 
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/  UP/DOV.N  INFERRED  RELATIVE  FAULT 

/DISPLACEMENT  ON  BLOCK 
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BOUNDARIES. 


A.  FAULTS  SHCV/N  ON  GEOLOGIC  CASE  MAP 


1  |  ALLUVIAL  MATERIAL 

|  |  ROCK  (ALL  PATTERNS) 


CONTOUR  INTERVAL  =  500FT. 
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DEPTH  CALCULATIONS  BASED  ON  DENSITY 
CONTRAST  OF  -0.5g  cm3 


GEOLOGIC  BASE  MAP:  STEWART  0  CARLSON  (1978) 
&  SELECTED  VERIFICATION  SEISMIC  REFRACTION  RESULTS 
BORINGS  FROM  LITERATURE 


DEPTH  TO  ROCK  INTERPRETED  FROM 
GRAVITY  DATA 
BIG  SMOKY  VALLEY,  NEVADA 
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